% 创建时间：2025-06-21 23:49:03
% 作者：兪馳
clc;
clear;

A=[-7,2,0,-1; 
    1,-4,2,1; 
    2,-1,-6,-1; 
    -1,-1,0,-4]; % 矩阵输入

syms t, A1=(expm(A*1j*t)-expm(-A*1j*t))/(2*1j); 
A1=simplify(A1) % 正弦函数
A2=(expm(A*1j*t)+...
    expm(-A*1j*t))/2; 

A2=simplify(A2) % 利用 Euler 公式求余弦函数

result1=simplify(rewrite(A1,'sin'));
result2=simplify(rewrite(A2,'sin'))
disp('---------------------------------------')
latex(result1)
disp('---------------------------------------')
latex(result2)